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We formulate empirically the rapidity density distribution of produced particles in multiple par- 
ticle production. The assumed mechanism is that the produced particles are emitted isotropically 
from several emitting centers, located on the rapidity axis. The formula includes five adjustable 
parameters, which are to be determined by the experimental data of (pseudo-)rapidity density dis- 
tributions and transverse momentum distributions at various energies. It is a distinguished difference 
of the present rapidity density distribution from those of other models that the particle production is 
suppressed strongly in the forward region. We discuss multiplicity and inelasticity at high energies, 
the pseudo-rapidity density distribution at LHC energy and some speculations, based on the present 
formulation. 



I. INTRODUCTION 

We formulate the rapidity density distribution of 
produced particles phenomenologically and analyti- 
cally on the basis of simple assumptions. It may 
clarify what kinds of mechanism are necessary essen- 
tially to describe multiple particle production and pro- 
vide a model which can be extrapolated with more 
confidence into higher energies. 



II. ASSUMPTIONS AND FORMULAE 

l. Assumed mechanism of multiple particle 
production 



Distribution of emitting centers 




FIG. 1: Produced particles are emitted isotropically from 
several emitting centers, located on the rapidity axis y* 
in CMS (the center of mass system). The distribution of 
emitting centers is a simple flat one of eq.(2) (the chain 
line). The Fermi distribution of eq.(3) (the chain-dot line) 
is shown together for comparison. y%a ax = \n(y/s/M) (y/s 
: the energy in CMS, M : nucleon mass) and yo = y^ax — 
lna2 (<i2 : an adjustable parameter). 



We assume the following for the mechanism of mul- 
5^ ' tiple particle production. 

(1) Produced particles are emitted isotropically from 
several emitting centers, which are distributed on the 
rapidity axis. (Fig. 1) 

(2) Produced particles are the newly produced ones 
excluding the surviving particle. We assume that 
all produced particles are pions (mass m) tenta- 
tively. (An effect of kaons among the produced par- 
ticles is not large for the (pseudo-) rapidity density 
distribution. [l[) 

(3) The normalized energy distribution of produced 
particles in the rest frames of respective emitting cen- 
ters is 



f(p)dp 



P" 



2(1 + r) 



1 
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-P/Pi 



i4 



~P/P2 



dp (1) 



where the values of the parameters p\ and pi {p\ < P2) 
are determined so as to reproduce the p T (transverse 



momentum) distributions at y/s = 63, 546 and 1800 
GeV. That is, the first and the second term in the en- 
ergy distribution correspond to the p T distributions in 
the low and high p T region, respectively. (We showed 
in Ref. [l| that the p T distribution cannot be described 
by a single exponential function of the energy distri- 
bution.) The parameter r is the relative contribution 
of the second term. (The term p 2 is necessary to re- 
produce the p T distribution in the vicinity of p T = 0.) 
(4) The distribution of the emitting centers in CMS 
(the center of mas system) is a "simple flat" one; 



g(y'W = 



yo 



(0 < y' < yo) 



(2) 



where y = y* max - lna 2 and y* max = lnjS/i/M) (y/s : 
the CMS energy, M : nucleon mass) . |20| (see Fig. 1) 
The parameter ai describes shrinkage of the 
(pseudo-)rapidity density distribution in the forward 
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region, and the value of it is determined so as to re- 
produce (pseudo-)rapidity density distributions in the 
region-3 2lj at various energies. Note that the distri- 
bution is normalized to one in the forward hemisphere. 

We examined another type of the emitting center 
distribution, 



g(y'W = 



;dy' 



(3) 



l + eHy'~y°) 

with b = 2.0 and c = 6/ln(l + e by °), which will be 
called "Fermi distribution", (see Fig. 1) 



where 



V* = In 



p T (& - e~ v ) + y/pU e ' f ~ e ~''*) 2 + V 



2//. 



dx*dp T y/(x*) 2 + (2Ai/v^) 2 dy*dp T 



where 



y* = In 



^/(x*) 2 + (2 M /Vi) 2 + x* 



(6) 



B. Rapidity density distribution 

The normalized energy-angular distribution of a 
produced particles in the rest frame of an emitting 
center is 

f(p)dp^d(cos9) 

where the variable is the zenith angle of the pro- 
duced particle in the rest frames of respective emit- 
ting centers. By the variable transformation it turns 
to 

-^dydp T 

where the variable y is the rapidity in the rest frames 
of respective emitting centers. Since the rapidity of a 
produced particle in CMS is y* = y' + y, the rapidity 
density distribution of charged produced particles in 
CMS is 



dy*dp T 



aiVo 



P -^f{p)g{y'W 
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y'), p = \/E 2 — m 2 and 
2 The parameter a\ is related to the 



where E = fx cosh(y 
fl = y/p 2 +m 
(pseudo-) rapidity density at y* = 0, and the value of 
it is determined so as to reproduce the data at various 
energies. Note that the proportional coefficient is not 
a\ but a±yo in eq.(4). 



D. p T distribution 

The p T distribution (at 6* — 90°) in terms of the 
invariant cross section is 



E 



d 3 p 



0~inel ( d 2 N ch 



>=90° 



2Ttp T \dy*dp 7 



(7) 



The suffix 8 = 90° will be omitted hereafter. 

The local p T average at the rapidity y* is defined as 



< Pt >y* : 



d 2 N ch 
Vt-t—, — dp 7 

o dy*dp T 



d 2 N ch 
dy*dp T 



dp-z 



(8) 

The local p T average at the pseudo-rapidity 
n* , < p T >r)*, is defined similarly, which is not the 
same as eq.(8) in general. 

Note that the p T average, obtained by the experi- 
ments, is either < p T > y *=o or < p T > v * = q in most of 
the experiments. 



E. Multiplicity and inelasticity 

Since the distribution of the emitting centers is nor- 
malized in the forward hemisphere, charged multiplic- 
ity m c h is twice that of the normalization coefficient 
in eq.(4), 



m ch = 2ai?/o = 2ai 



ha( — 1-lnoa 



(9) 



The total inelasticity in CMS is 



Pseudo-rapidity density distribution and 
^-distribution 



K=- 
2./o 



*• ( ^ I dx* 
V dx* 



(10) 



The pseudo-rapidity density distribution and x- 
distribution (defined as x* = 2p*/y / s, p* : longitu- 
dinal momentum of produced particle in CMS) are 
obtained by variable transformations from eq.(4). 



d 2 N, 



ch 



d 2 N, 



ch 



drj*dp T ^Jp 2 (e^* + e~^) 2 + Am 2 dy*dp 7 



(5) 



III. VALUES OF THE PARAMETERS 

A. Values of the parameters pi, pi and r 

The values of the parameters in the energy dis- 
tribution of cq.(l) arc determined by fitting the p T 
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TABLE I: Values of the parameters for the p T distributions 



(GeV) 


Pi or p2 
(GeV/c) 


(dN ch 
cal. 


/drj*) v *=o 
exp. 


Ol 


(7 incl 

(mb) 


r 




Ref. 


63 


pi = U.104 

p 2 = 0.632 


0.843ai 


1.88 ±0.08 


2.23 


36.0 


2.0 x 10" 


-4 


[2] 


546 


pi =0.1175 
p 2 = 0.895 


0.865ai 


2.79 ±0.08 


3.23 


49.0 


2.0 x 10" 


-3 


[2] 


1800 


pi = 0.206 
p 2 = 0.938 


0.875ai 


3.95* 


4.51 


56.0 


2.0 x 10" 


-3 


[3] 


1.4 x 10 4 


pi = 0.253 
p 2 = 0.632 






6.54 


73.0 1 " 


6.5 x 10" 


-2 





* : the value at rj* = 0.12, * : Model 2 in Ref. [7] 




P T (GeV/c) 



FIG. 2: The p T distributions at y/s = 63 GeV (the lower 
curve) and 546 GeV (the upper). Experimental data are 
compiled by UA1 Collaboration. [2] 

dots ((h+±h")/2) : UA1 Collaboration for a half of 
charged hadrons at y/s — 546 GeV, squares ((h + ±h~)/2), 
diamonds (tt°), triangles (ty°) and circles ((h + ±h~)/2) : 
ISR data at y/s = 63 GeV. Note that value of the curve 
in the ordinate is a half of eq.(7) with the values of the 
parameters in Table 1. 



distribution of eq.(7) to those of the experiments at 
y/s = 63, 546 and 1800 GeV, HQ which are shown in 
Figs. 2 and 3. 

Values of the parameters determined are tabulated 
in Table 1. The values of the parameters p\ and pi 
are determined to reproduce the slopes of the p T dis- 



tributions in the p T regions of 0.2 — 1.0 and 5.0 — 10.0 
(GeV/c), respectively. The value of the parameter 
ai in Table 1 is determined by equating the calcu- 
lated pseudo-rapidity density at 77* = to that of the 
experiment. [<j (In the calculation we assume the value 
of the parameter = 1.0, which affects slightly the 
pseudo-rapidity density at n* = 0.) One can see in the 
figures that the p T distributions are described well in 
the small and large p T regions, but not in the middle 
p T region. 

Fig. 4 shows the energy dependences of the param- 
eters pi, pi and r in Table 1. Assuming the power 
dependences of the energy, empirical formulae of the 
parameters are; 

pi = 0.0895(Vs) - 109 

P2 = 0.381U/I) - 130 (11) 

r = 2.38 x 10- 6 (Vs) 107 

where the energy y/s is in GeV. 

We expected that the parameters p\ and P2 are 
energy-independent, which is not correct. The value 
of the parameter r exceeds 1.0 for the energy y/s > 
1.80 x 10 5 GeV, since the exponent is as large as 1.07. 
We can append, however, that the cross section of 
mini-jets (E T > 5 GeV, \n*\ < 1.5) increases simi- 
larly in the energy region y/s — 200 — 900 GeV. [5| We 
should note that the exponent may become smaller 
than 1.07 at high energies, since the channel of mini- 
jets has opened just in the present energy region. 

Fig. 5 shows the p T average at n* — 0, < p T >^*=0j 
and the energy y/s, based on the energy dependences 
of the parameters in eq.(ll). It is of no wonder that 
experimental data are consistent with that at 77* = 
(but not at y* = 0), since the condition to sample 
the events refers the pseudo-rapidity in the data con- 
cerned. 

Fig. 6 shows the local p T average at the rapidity 
y* for various values of the parameter 02, together 
with the data in the region-3. The data are from 
UA7 Collaboration at y/s = 630 GeV.Q The data 
are described well by the curve of 02 = 5.0. It is im- 
portant to note that the rapidity density distribution 
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FIG. 3: The p T distributions at 1800 GeV (the lower 
curve) and 14 TeV (the upper). Experimental data are 
for a half of charged hadrons, (h + +h~)/2, from CDF Col- 
laboration. [3] Note that value of the curve in the ordinate 
is a half of eq.(7) with the values of the parameters in 
Table 1. 




Energy vTs" (GeV) 

FIG. 5: The p T average and the energy y/s. The data are 
those compiled by UA5 Collaboration[2] and by UA1 Col- 
laboration [8] and those from CDF Collaboration. [3] The 
curves are < p T > y *=o (denoted as "a") and < p T >,,*=o 
("b"), which are defined in 2.4. 




Rapidity y* 

FIG. 6: The local p T average of eq.(8) at y/s — 630 
GeV. The curves are for the values of the parameter 
ai — 1.0, 2.0, 3.0, 5.0 and 7.0 (attached to the curves). 
Data are for 7r°'s by UA7 Collaboration. [6] The data are 
described well by the curve of a% = 5.0. 



10 



0.1 




10 100 1000 5000 



Is (GeV) 

FIG. 4: Energy dependences of the parameters p\ (GeV/c) 
(open circles), p2 (GeV/c) (full circles) and r (squares). 
Note that the parameter r is multiplied by a factor 10 3 . 
The solid lines are eq.(ll), the best-fit ones to the data 
points, assuming the power dependences of the energy y/s. 



in the region-3, obtained by the same collaboration, 
is described by the same value of the parameter a-z 
simultaneously, (see 3.2) 



B. Values of the parameters ai and a 2 

The values of the parameters a\ and are deter- 
mined by fitting the (pseudo-)rapidity density distri- 
bution, eqs.(4) or (5), to those of experiments at var- 
ious energies. Note that the surviving particle is in- 
cluded among the observed particles in most of the 
experimental data. The density of the surviving par- 
ticle, however, occupies a small part of the (pseudo-) 
rapidity density in the region-1 and -2 where the data 
exist, (see, for example, Fig. 13) 

Sources of the experimental data of the (pseudo-) 
rapidity density distribution are tabulated in Table 2. 
The data in the table concern only those that are nec- 
essary for the present analysis. For example, EHS- 
NA22 Collaboration observed K + p and ir + p collisions, 
too. Necessary comments to respective sets of data are 
found in Rcf. 1]. 

Figs. 7, 8 and 9 are the examples to show how well 
the experimental data are described by the curves of 
the present formulation. The curves arc calculated 
for several assumed values of the parameters a\ and 
et2, i.e. a± = 1.0 and ai — 1.0, 2.0, Hence, in 
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TABLE II: Data of (pseudo-)rapidity density distributions and the values of the parameters ai and a-z 







Energy 




Observed Observed 




Observed 


Parameter 




Collab. 


Site 


yfs (GeV) Collision particles* 


range 




quantity 


ai 


a 2 


Ref. 


EHS-NA22 CERN SPS 


22.4 


PP 


7T , C + 


\V*\ < 6-0 




da/dn* 


1.8 


1.5 ±0.5 


[9] 


Phobos 


BNL RHIC 


200 


PP 


C± 


\V*\ < 5-4 




dN/dn* 


2.6 


4± 1 


[10] 


UA5 


CERN SPS 


53 


PP 


C ± 


\V*\ < 3-5 




dN/dn* 


2.1 


3± 1 


W 

I J 






200 






\V*\ < 4 - 6 






2.6 


4± 1 








546 






< 4-6 






3.1 


8.5 ± 1.5 








900 






\t}*\ < 4.6 






3.5 


10 ±1 




UA7 


CERN SPS 


630 


PP 


7 


y* = 5.0- 


6.6 


da^a/dy* 


(3.5) 


5± 1 


[6] 


P238 


CERN SPS 


630 


PP 


c± 


n* = 1.5- 


5.5 


dN/dn* 


3.5 


6± 1 


[11] 


CDF 


FNAL Tevatron 


630 


PP 


c± 


\V*\ < 3-5 




dN/dn* 


3.6 




[12] 






1800 












4.4 







The letter c stands for charged particles. 



the figures, the calculated curves are shifted upwards 
to fit to the data in the region- 1, and then the curve, 
which fits best to the data in the region-2 and/or -3, is 
selected. This procedure determines the values of the 
parameters a\ and ai, which are tabulated in Table 2. 
(The value of the parameter ai cannot be determined 
for the data by CDF Collaboration, since the data exit 
only in the region- 1.) We list some comments below 
to the figures. 




Pseudo-rapidity r[* 

FIG. 7: Pseudo-rapidity density distribution at *J~s = 22.4 
GeV. The data are the sum of 7v~ (negative pions) and c + 
(positive particles) by EHS-NA22 Collaboration. [9] Dotted 
curves are eq.(5) with the values of the parameters a\ = 
1.0 and <i2 = 1.0, 2.0 and 3.0 (attached to the curve). 
The solid lines are the best-fit ones to the experimental 
data with the values of the parameters a\ = 1.8 and ai = 
1.0 — 2.0. Agreement between the curves and the data is 
poor in the region-2. 

(1) In Fig. 7 the agreement between the best-fit curves 
and the data is good in the region-1 and -3, but is 
poor in the region-2. This situation is the same for al- 
most all sets of data in Table 2. In Ref.[l| we showed 
that the Fermi distribution for the emitting centers, 
eq.(3), brings better agreement in the region-2, but 
that it does not describe the data of the local p T aver- 
age in the region-3, obtained by UA7 Collaboration. 
One can see, however, in 3.1 that the simple flat dis- 



tribution for the emitting centers describe the local p T 
average naturally, which is the reason why we adopt 
it in the present paper. 



0.01 







UA5(546) j 
a,= 7-10 • 








10 7 5 


\ 1 [ 



Pseudo-rapidity r|* 



FIG. 8: Pseudo-rapidity density distribution at ^/s = 546 
GeV. Data are by UA5 Collaboration. [4] Dotted curves are 
eq.(5) with the values of the parameters a\ = 1.0 and <i2 = 
1.0, 5.0, 7.0 and 10.0 (attached to the curve). The solid 
lines are the best-fit ones to the experimental data with the 
values of the parameters oi = 3.1 and a-2 = 7.0—10.0. The 
chain line is for the values of the parameters a\ = 3.1 and 
<i2 = 5.0, which corresponds approximately to the assumed 
case of the inelasticity K = 0.5. (see eq.(12)) Note that 
the chain curve does not describe the data clearly. 

(2) Fig. 8 shows the pseudo-rapidity density distri- 
butions at ^/s = 546 GeV, together with the experi- 
mental data by UA5 Collaboration. [4| The curves of 
ai =3.1 and ai — 7 and 10 describe the data reason- 
ably. 

(3) In Fig. 9 the value of the parameter a\ = 3.5 for 
the data by UA7 Collaboration Q is the assumed one, 
since data exist only in the region-3. It is important to 
note that the rapidity density distribution is described 
well by the value of the parameter ai — 5.0 which 
describes the local p T average in Fig. 6. 

Fig. 10 shows the energy dependences of the param- 
eters ai and ai in Table 2. We list some comments to 
Fig. 10. 
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FIG. 9: Rapidity density distribution at yfs = 630 GeV. 
Data of 7r°'s, by UA7 Collaboration [6], are multiplied by 
a factor 2.0. Dotted lines are those for the values of the 
parameters ai = 1.0 and ai = 1.0, 2.0, 3.0, 5.0 and 7.0 
(attached to the curves). The solid line is the best-fit one 
with the values of the parameters a\ = 3.5 (assumed) and 
a 2 = 5.0. 



FIG. 10: The values of the parameters a 1 (open 
marks) and ai (full marks) in Table 2 and the energy 
yfs. Marks: triangles (EHS-NA22 Collaboration), cir- 
cles (UA5), squares (Phobos), pentagons (P238), dia- 
monds (UA7), stars (CDF). Solid lines are the best-fit 
ones to the data points, eq.(13) in the text. Chain lines 
are those in Ref. [1]. The chain-dot line is eq.(12) in the 
text, for which the total inelasticity is K — 0.5. 



(1) Since the (pseudo-)rapidity density distribution is 
almost flat in the region- 1, it is easy to fit the data to 
the calculated curve. Consequently the values of the 
parameter a\ are determined reliably to reproduce the 
data. 

(2) In the accelerator experiments it is not easy to 
obtain the data in the region-2 and -3 due to the ex- 
perimental conditions. Consequently the data in the 
concerned regions are limited or missed in most of the 
experiments in Table 2, and the data in the region- 
2, if they exist, are often not consistent one another 
by the experiments even at the same incident energy. 
Consequently the values of the parameter a 2 are de- 
termined with large errors and are distributed widely. 

(3) In order to consider the assumed case of the inelas- 
ticity K — 0.5, we examine the energy dependence of 
the parameter 

a 2 = 0.718( % /i) ' 320 (12) 

As can be seen in Fig. 12, eq.(12) brings the inelastic- 
ity K ~ 0.5, taking the value of the parameter a\ as 
it is, one of eq.(13). In Fig. 10 the chain-dot line of 
eq.(12) is almost consistent with points except those at 
yfa = 546 and 900 GeV, both of which are from UA5 
Collaboration. (Fig. 8 shows that the pseudo-rapidity 
density distribution at yfs — 546 GeV by UA5 Collab- 
oration is not described by the curve of the assumed 
case of K = 0.5 clearly.) It is not evident, however, 
that both data are biased in the region-2 seriously. 
Since it is our strategy in the present paper to formu- 
late multiple particle production phenomenologically 
avoiding a priori assumptions as much as possible, we 
determine the energy dependence of the parameter a 2 
by the least square method including both data. 

Assuming the power law for the energy dependences 
of the parameters a% and a 2 , we obtain (yfs in GeV) 

ax = 0.915(Vi) - 206 n ^ 
a 2 = 0.289(v/i) - 501 1 ' 



which are shown in Fig. 10. 



IV. DISCUSSIONS 

(i) Energy dependence of the multiplicity 

Fig. 11 shows the energy dependence of the charged 
multiplicity, eq.(9), by the present formulation, to- 
gether with experimental data. [|[ It describes the data 
well at high energies of yfs > 50 GeV. The data at 
lower energies are described better by the curve of the 
assumed case of the inelasticity K = 0.5. (see the text 
below for the value of the inelasticity at low energies) 

LHC l(f eV 




10 l(f 10' io 4 10' 10° 

Energy >Ts~ (GeV) 



FIG. 11: Energy dependence of the charged multiplicity 
m c h, defined by eq.(9). The solid line is for the values of 
the parameters in eq.(13). Experimental data are compiled 
by UA5 Collaboration. [4] The chain line is for the assumed 
case of eq.(12), for which the inelasticity remains constant 
(-0.5). 



(ii) Energy dependence of the inelasticity 

Fig. 12 shows the energy dependence of the total 
inelasticity, eq.(10), by the present formulation. The 
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together in Fig. 13. We assume that the charge ex- 
change probability of the incident proton into neu- 
trons is 0.5, details of which are described in Ref.[l[. 



L:niT!jy 



FIG. 12: Energy dependence of the total inelasticity, de- 
fined by eq.(10), by the present formulation, (the solid 
line) The inelasticity decreases with the energy in the en- 
ergy region ^fs < 2 x 10 4 GeV, and then increases. 

The shaded area is the predictions by some models of 
multiple particle production. [13] Note that the predicted 
quantities are not exact inelasticities but (1 — r/)'s where 
the parameter n' is the energy fraction of the highest en- 
ergy baryon among the produced particles. ("Produced 
particles" here include the surviving particle, which is dif- 
ferent from the definition in the present paper.) 

The dotted line is ai/ai oc (\fs)~ ' 295 , which is steeper 
than the inelasticity due to the energy dependence of the 
parameters pi and p2 in the energy distribution. The chain 
line is for the assumed case of eq.(12), for which the in- 
elasticity remains almost constant (~ 0.5). 



inelasticity decreases with the energy in the energy 
region s/s < 2 x 10 4 GeV, due to the rapid increase 
of the parameter a-i, compared with the parameter 
a\. (The inelasticity is proportional to a,\jai approx- 
imately.) After that it increases, due to the exchange 
of the dominant parameter in the energy distribution 
from pi to p 2 (pi < P2) through rapid increase of the 
parameter r. (Experimental data of the inelasticity 
is K ~ 0.5 at low energies around y/s — 10 GeV by 
bubble chamber experiments.) 

(hi) Pseudo-rapidity density distribution at LUC en- 
ergy 

Fig. 13 shows the pseudo-rapidity density distri- 
bution of the produced particles at = 14 TeV 
(LHC energy) by the present formulation. Compared 
with the predictions by other models, the shrinkage of 
the forward region is distinguished, which results in a 
small inelasticity of the present formulation. Predic- 
tions by the models are distributed widely, and even 
a single data of the pseudo-rapidity density at r\* = 
can discriminate some models as improbable if they 
cannot modify their predictions. [22| 

According to CMS Collaboration the pseudo- 
rapidity density at if = at \fs = 7 TeV is 
dN ch /dri* = 5.78 ± O.Ol(stat) ± 0.23(syst).[H| It in- 
dicates that the increase of the parameter a\ is more 
rapid than the one of eq.(13) at s/s = 7 TeV, which 
will be discussed elsewhere. 

The distribution of the surviving particle is shown 



Q 




= Uiltf (GeV) 








P : x \\ 















































Pseudo-rapidity T|* 

FIG. 13: Pseudo-rapidity density distribution of produced 
particles (the solid line) and the surviving particle (the 
shaded area) at \/s = 14 TeV (LHC energy). The thin 
lines are the predictions of the pseudo-rapidity density by 
some models. [14, 15] 

Q : QGSJET01, S : SIBYLL2.1, D : DPMJET2.55, 
P : PYTHIA, n : neXus2.1 



(iv) Speculations related to the inelasticity to decrease 
and to increase 

According to the present formulation the inelastic- 
ity is decreasing in the energy region s/s < 2 x 10 4 GeV 
(E < 2 x 10 17 eV) and then increasing in \/s > 5 x 10 4 
GeV (E Q > 10 18 eV). This structure of the inelastic- 
ity is caused by the rapid increase of the parameter 
<22, compared with that of a\, which means that par- 
ticle production is suppressed strongly in the forward 
region. 

The inelasticity is related to the attenuation mean 
free path of cosmic rays, which is a dominant factor to 
govern the cosmic-ray propagation in the atmosphere. 
Small (Large) inelasticity makes the attenuation mean 
free path long (short). Consequently the development, 
rise and fall, of the air showers becomes slow (rapid), 
and the air shower size at the maximum development 
is small (large) since the total track length of the air 
shower particles is conserved. Hence there is a possi- 
bility that following problems may be cleared by the 
decreasing and increasing inelasticity. 

(1) Intensity of the primary cosmic rays. 

There is a discrepancy between the primary cosmic- 
ray intensities by balloon experiments of direct obser- 
vation and by air shower experiments of indirect ob- 
servation, the former being lower than the latter by a 
factor 2.[l7| 

(2) < X max > and RMS(X max ) by Auger Collabora- 
tion 

According to Pierre Auger Collaboration to observe 
highest energy air showers, both the depth of shower 
maximum, < X max >, and the dispersion of the max- 
imum depth, RMS(X max ), are reaching the expected 
lines of the iron primaries from those of the proton 
primaries in the region E = 10 18 — 4 x 10 19 eV[18| 
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Detailed and quantitative discussions will be made elsewhere. 
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